Deconfined criticality, runaway flow in the two- component scalar 
electrodynamics and weak first-order superfiuid-solid transitions 

A. B. Kuklov\ N. V. Prokof'ev^'^.^, B. V. Svistunov^'^, and M. Troyer^ 
^Department of Engineering Science and Physics, 
The College of Staten Island, CUNY, Staten Island, NY 10314 
"^Department of Physics, University of Massachusetts, Amherst, MA 01003 
^ Dipartimento di Fisica, Universitd di Trento and 
CNR-INFM EEC Center, 1-38050 Povo, Trento, Italy 
'^Russian Research Center "Kurchatov Institute", 123182, Moscow 
^ Theoretische Physik, ETH Ziirich, CH-8093 Zurich, Switzerland 

Abstract 

We perform a comparative Monte Carlo study of the easy-plane deconfined critical point (DCP) 
action and its short-range counterpart to reveal close similarities between the two models for 
intermediate and strong coupling regimes. For weak coupling, the structure of the phase diagram 
depends on the interaction range: while the short-range model features a tricritical point and 
a continuous U(l)xU(l) transition, the long-range DCP action is characterized by the runaway 
renormalization flow of coupling into a first (I) order phase transition. We develop a "numerical 
flowgram" method for high precision studies of the runaway effect, weakly I-order transitions, and 
polycritical points. We prove that the easy-plane DCP action is the field theory of a weakly I-order 
phase transition between the valence bond solid and the easy-plane antiferromagnet (or superfluid, 
in particle language) for any value of the weak coupling strength. Our analysis also solves the long 
standing problem of what is the ultimate fate of the runaway flow to strong coupling in the theory 
of scalar electrodynamics in three dimensions with U(l)xU(l) symmetry of quartic interactions. 

PACS numbers: 
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INTRODUCTION 



Critical properties of systems described by several complex fields coupled to the gauge 
field have a long history of studies and are relevant for numerous problems in physics which 
include normal-superfluid transitions in multi-component neutral or charged liquids (see, e^. 




I), superfiuid-valence-bond solid (SF-VBS) transitions in lattice models 
the Higgs mechanism in particle physics 8], etc. Recently, the authors of Refs. 
argued that the SF-VBS transition in a (2 + l)-dimensional system is an example of a 
qualitatively new type of quantum criticality ( "decondined criticality" ) that does not fit the 
Ginzburg-Landau- Wilson (GLW) paradigm. As suggested in reference an experimental 
realization of the SF-VBS transition is possible in the system of ultracold atoms trapped in 
optical lattice. If true, the significance of the claim that DCP opens a new era in the theory 
of phase transitions is hard to overestimate. The study of the SF-VBS transition within 
the DCP theory framework, and the closely related problem of the runaway flow to strong 
coupling in the theory of scalar quantum electrodynamics in three dimensions (3D) P, 
is the main focus of this work. 

We also discuss a short-range analog of the DCP action, namely, the model describing 
lattice phases of the two-component Bose-Hubbard (BH2) system with large integer filling 
factors [l]]. The only difference between the DCP and BH2 models lies in the interaction 
range. The DCP action is characterized by a long-range Coulomb type interaction, which 
is replaced with the contact (5-functional form in the BH2 analog. Despite this difference, 
the gross features of the phase diagram are surprisingly similar in both models, especially 
in the regimes of intermediate and strong coupling between the two field components. Even 
quantitatively, the phase diagrams are nearly indistinguishable, see Fig. |21 Both feature a 
bicritical point (BP) above which there are three phases: a superfiuid (SF), a paired phase 
and an insulator. The superfiuid is characterized by the order parameters for two complex 
scalar fields, {ipi) ^ 0, {1IJ2) 7^ 0. In the paired phase the order parameter is {ipl4'2) 7^ 
while (ipi) = {1P2) = 0. In the context of the DCP theory, the paired phase is associated 
with the valence-bond supersolid (SFS); within the BH2 model it corresponds to the 
super-counter-fluid (SCF) phase describing pairing between particles of one component and 



holes of another 



12| . Finally, the insulating phase represents the VBS state in the DCP 



theory and the Mott insulator (MI) state in the BH2 model. In this phase all off-diagonal 
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order parameters are zero. 

The crucial difference between the two models is seen below the bicritical point BP, where 
direct transitions from insulator to SF take place. In the DCP action this line describes 
the SF-VBS ..a„*™a«o„. An i.p.,* a„.u,„p«o„ ™ade in Refs. fl B Q was that as 
the interaction constant between the field components, g, weakens down to g = (where 
DCP and BH2 actions become identical and describe two decoupled XY-models), there 
exists a lower tricritical point (TP) giTP > below which the VBS-SF transition becomes 
continuous. Since lower TP does exist in the short-range BH2 model the expectation 
was that both models have essentially identical phase diagrams, with only one important 
difference — in the short-range case the continuous transition below lower TP accounts for the 
U(l)xU(l) universality, while the continuous transition in the DCP action for < (7 < Q ltp 
would be in the new "deconfined" universality class. However, as we suggested earlier [l3| 
and now clearly demonstrate below, no such lower TP exists in the DCP action and the line 
of I-order SF-VBS transitions extends all the way to g ^ 0. 

The tool we introduce for studying subtle features of the phase diagram is the numerical 
flowgram method. It is based on scaling properties of quantities similar to Binder cummu- 
lants Ql and turns out to be much more adequate for analyzing weakly I-order transitions 
than the conventional method based on the bimodal shape of the energy/action distribution 
and hysteresis loops. The flowgram technique can also be used for systematic identification 
of polycritical points. To the best of our knowledge no other method can achieve this goal 
with comparable accuracy. First, we construct the flowgram for the BH2 model and demon- 
strate its efficiency for those parts of the phase diagram which can be reliably identified by 
more conventional methods. Next, we apply the new method to the most controversial re- 
gion of the BH2 phase diagram between BP and upper TP, which we failed to resolve in the 
previous study 11]. The existence of upper TP proves that the mean field does reproduce 
the topology of the phase diagram for the short-range model correctly, with the reservation 
that fluctuation-induced effects significantly shorten the interval between the bicritical and 
the upper tricritical points. 

Finally, we apply the flowgram method to the self-dual version of the DCP action P|. 
We have identified the bicritical and the upper tricritical points, but no lower TP has been 
found. We directly observe the runaway flow to strong coupling and prove that its ultimate 
fate is a I-order transition. More specifically, for small g, we observe data collapse with 
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respect to the coupling strength rescaled by the system size L. This collapse indicates the 
existence of an effective long-range coupling ges{L) and proves that it is never weak in the 
L ^ oo limit: no matter how small is g at short scales it always gets renormalized to Qcs ~ 1 
where the I-order transition takes place. Perfect data collapse continuously relating small g 
to one for which I-order transition is observed at accessible sizes provides crucial evidence 
against lower TP and in support of the claim that the direct VBS-SF transition is always 
weakly I-order in accordance with the Ginzburg-Landau- Wilson paradigm [iS^ . 

On one hand, our results represent the first solution of the long-standing runaway problem 
for the two-component field-theory in question. On the other hand, they explain the outcome 
of recent numerical simulations which report either difficulties with finite-size scaling of the 
data |15] or weakly I-order transitions in microscopic models of the SF-VBS transition. 



DCP ACTION AND ITS SHORT-RANGE COUNTERPART 

There are many equivalent formulations of the action describing multi-component scalar 
fields coupled by the Abelian gauge field . In what follows we will concentrate on the 

case of two identical non-convertible components and will employ the integer-current lattice 
representation which can be viewed as either a high-temperature expansion for XY models in 
three dimensions, or as a path-integral (world-line) representation of the interacting quantum 
system in discrete imaginary time in {d + 1) = 3 dimensions. The DCP action is explicitly 
symmetric with respect to exchanging the components. After the gauge field is integrated 
out, it reads 

S? = + Ujl + gQr-r' (jlr + j2r) " (jlr' + hv) , (l) 

with U > 0, g > 0. The lattice Fourier transform of the long-range interaction po- 
tential Qr-r' is given by Qq = 1/ sin^(g^j/2), which implies the asymptotic form 
Q{r —>■ oo) ~ 1/r. The summation over sites of the simple cubic space-time lattice with 
periodic boundary conditions, r, is assumed. Integer currents j^r = {jar)fi with a = 1,2 and 
fi = x,y,T are defined on the lattice bonds and are subject to the zero-divergence, Vj = 0, 
constraint, i.e. the configuration space of j-currents is that of closed oriented loops, see 
Fig. H In terms of particle world lines, currents in the time direction represent occupa- 
tion number fiuctuations away from the commensurate filling, and currents in the spatial 
directions represent hopping events. In the discussion of the SF-VBS transition, ji and j2 
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FIG. 1: (Color online). A typical configuration of currents for the two-component j-current model. 
The configuration space consists of arbitrary number of closed oriented loops which are allowed to 
overlap and intersect. The current ja=i,2 is defined as an algebraic sum of all same-type arrows on 
the corresponding bond. 



currents in the DCP action represent world lines of spinons which are VBS vortices carrying 
fractional particle charge ±1/2 0,0]. The insulating VBS state is characterized by small 
current loops; the corresponding spinon order parameter, tpa, = 1? 2, is zero in this phase. 
When spinon world lines proliferate and grow macroscopically large, the system enters the 
superfluid state. The paired phase (SFS) occurs at large g, when only loops in the ji — j2 
channel grow macroscopically large. Since bound pairs of spinons with opposite vorticity 
make particle and hole excitations, the SFS state is a superfluid with the solid VBS order, 
i.e. a supersohd. 

As discussed above, many features of the DCP action are remarkably similar to those of a 
more conventional two-component Bose- Hubbard model with the short-range interaction 




FIG. 2: (Color online). Phase diagrams of the long-range (left panel), Eq. and short-range 
(right panel), Eq. actions. Bold solid lines indicate the I-order phase transition. Error bars 
are shown but are typically much smaller than symbol sizes. 



potential 

Si'^ = Uil + U}l + V ihr+hr? , (2) 

with U > 0, V > 0. Phase diagrams for the long- and short-range actions are shown in 
Fig. 121 Different phases are identified in terms of loop sizes for ji and j2 currents similar 
to the discussion given above: the Mott insulator (MI) is characterized by small ji- and 
j2-loops, in the two-component SF there are macroscopically large ji- and j2-loops, and in 
the SCF state only ji — j2 loops are macroscopically large. 

Previous Monte Carlo data for the BH2 model Q did not agree with the mean field 
theory wh.ch describes the system in terms of three fields - for each cornponerrt and 
$ for the paired field Omitting gradient terms the mean-field action is given by: 

smf = Ih\h' + + r^mm'] + limH' + + uM'] 

- 7($A>2 + c.c.) , (3) 
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where r, Vm, > 0, f/^ > 0, 7 are some effective coefficients which depend on the bare 
couphng parameters U, V in the action (0)- Trivial minimization of Q with respect to the 
order parameter fields results in the phase diagram which can be formulated in terms of 
two dimensionless variables ~ r and ~ r^, with all other coefficients rescaled to become 



unity A sketch of the diagram is presented in Fig. El It features a continuous SCF-MI 
transition, which terminates at BP point, as well as I-order MI-SF and SCF-SF transitions. 
The SCF-SF transition starts as a first-order line at BP and goes up to upper TP. Then, it 
continues as a second (II) order line for larger values of V. Similarly, the MI-SF transition 
is of I-order type between BP and lower TP. Then, it continues as a Il-order U{1) x U{1) 
transition line towards the = point where the components become decoupled. The 
length of the I-order line along the SCF-SF boundary is predicted to be about two times 
shorter than that along the MI-SF line, see Fig. El 

It is worth stressing that the I-order transition is induced by fluctuations of the paired 
field $. Indeed, one can approximately integrate $ out using gaussian approximation for 
this field and generate an effective interaction term ~ — (7^/|rm|)|^iP|'?/'2p- It destabilizes 
the system for sufficiently small \rm\ fill- Thus, the vicinity of BP must feature a I-order line 
(along MI-2SF and SCF-2SF). Correspondingly, there are two minimal values of {'y'^ /\rm\)m.in 
- one for > and the other for < - below which the transition becomes continuous 
along the respective boundaries. It is also important to realize that at the mean-field level 
there is no difference between the BH2 and DCP models. 

Monte Carlo simulations did find an extended I-order MI-SF boundary and lower 
TP, but failed to resolve conclusively BP and upper TP. Standard numerical tools based 
on studies of hysteresis and probability distributions for the action were not precise enough 
for this system (close to TP the I-order transition is very weak and hard to identify in 
finite-size simulations). In all other respects, the mean field theory and simulations agreed 
on the topology of the phase diagram. Another important observation is that simulations of 
the same model ^ in D = 3 + 1 dimensions [v}\ revealed an extended region of first-order 
transitions above BP. Since D = 4 is the upper critical dimension where fluctuations become 
essentially insignificant this result was expected. Thus, an intriguing possibility has emerged 
that the predicted first-order transition is destroyed by fluctuations in D = 2 + 1. 

In this work we apply a new flowgram method to perform a more refined analysis of the 
phase diagram in the vicinity of BP point and find that the SCF-SF boundary does contain 
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FIG. 3: (Color online). Sketch of the mean field phase diagram for the effective action (jSJ. The 
bold (red) solid line represents I-order transitions. Dotted lines are used for Il-order transitions. 
The top of the solid line corresponds to upper TP, and the bottom - to lower TP. All three phase 
boundaries meet at BP. 



a I-order region in D = 2 + 1. The distance between BP and upper TP turns out to be 
an order of magnitude smaller than predicted by the mean-field theory. We attribute this 
to the fluctuation induced suppression of the pairing effect so that higher V is required to 
produce the SCF phase. 

At this point one may ask why not to trust the mean-field theory on another prediction, 
namely, the existence of lower TP for the DCP action. The precise answer lies in how the 
effective coefficients in Q depend on U, g, and L. For the short-range model, the relation is 
rather straightforward, 7 ~ VV, leading to finite value of lower TP at some V = Vltp 7^ 0. 
Ostensibly, the long-range nature of interaction in (Q), which culminates in the runaway flow 
makes the relation between 7 and g system size dependent through ges = gL (in leading 
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order in g) and precludes simple perturbative reasoning for lower TP even for small g. In 
other words, pairing correlations are never weak at large scales in the DCP action, and it is 
possible that this model is always above lower TP of the mean field theory Q. 

FLOWGRAM METHOD 

In this Section we introduce the flowgram method and apply it to the study of the 
short-range model The most fundamental property of current loops responsible for the 
superfluid response of the system are winding numbers 



Integer Wa numbers count how many times jar currents wind around the system with periodic 
boundary conditions. The superfluid stiffness, A^, is directly proportional to (W^)/L. For 
any scale- invariant transition there should be ~ 0(1) loops of size i? in a volume R^, 
which implies that (W^) ~ 0(1) at criticality and directly leads to the Josephson relation 
Aa ^ l/L. Moreover, at the continuous transition point not only the (W^) value but also the 
whole distribution function F(W) is universal since it is governed by the system behavior 
at the largest scale. 

Given close similarities between the DCP and BH2 phase diagrams it is convenient to 
introduce the following terminology: the transition line between phases with zero and non- 
zero W_ = Wi — W2 winding number fluctuations is called a neutral line (described by 
proliferation of neutral j_ = ji — j2 loops), and the transition line between phases with 
zero and non-zero W_|_ = Wi + W2 winding number fluctuations is called a charged line 
(described by proliferation of charged = ji +j2 loops in addition to the loops of j_ which 
are already macroscopic). This terminology originates from the identification of j-currents 
as world line trajectories for charged spinon fields in the DCP action; only j_(.-currents are 
coupled by the long range Coulomb potential. Using this terminology, the SCF-SF (SFS-SF) 
transition is a charged line and the MI-SCF (VBS-SFS) transition is the neutral line; the two 
lines merge at BP. Notice that above BP the neutral winding number fluctuations (W^ ) on 
the charged line must exhibit linear growth with L because the paired phase is characterized 
by finite stiffness (see Fig. The same language will be used for the short-range action 
(j21) even though no Coulomb potential is involved. 




(4) 



r 
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Formally, in Eq. (^J only W_ numbers may have non-zero values in the thermodynamic 
limit. Indeed, charged winding numbers couple to Qq=o = = oo and thus are forbidden. 
However, charged loops of macroscopic size which simply do not wind around the system 
(but are similar to winding loops otherwise) are allowed and their contribution to the action 
can be estimated as Qci=i/L/L ~ L (one would get the same result for loops with non-zero 
winding numbers if the zero- momentum compotent of the interaction potential is "regular- 
ized": Qq=o <5q=27r/L ~ V-^^)- Slucc (5q=o = oo imposcs a global constraint on a single 
variable of the configuration space of the model it may not have any effect on the nature of 
the scale-invariant transition. For example, on approach to the critical point charged loops 
of larger and large size (but smaller then L — > oo ) will proliferate anyway. We remove the 
global constraint by setting (5q=o = Qq=27r/L for only one reason — to use non-zero winding 
number fiuctuations as an indication of the deconfinement transition. We expect (W^) ~ 1 
at the transition point. The advantage of looking at the W+ numbers is that they are 
sensitive only to the deconfinement transition to the SF phase and remain zero in both the 
VBS and SFS phases. 

The key elements of the fiowgram method are: 

• introduce a definition of the critical point for finite-size systems consistent with the 
thermodynamic limit and insensitive to the transition order. 



at the transition point, calculate quantities which are scale- invariant for the continu- 
ous phase transition in question, vanish in one of the phases and diverge in another; 
in our case such quantities are statistical averages of winding numbers squared. 



• study the fiow of R^{L) with system size L for given model parameters. As far as the 
thermodynamic critical points are concerned, the critical values of the corresponding 
parameters in the action are naturally obtained by extrapolating the sequence of finite- 
size critical points to L = cxo. 

Observing how R{L) scales with L at the transition point allows studying the order of a 
transition, since R{L oo) has to saturate at some universal constant for second-order 
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and to diverge for first-order transitions. Moreover, radical changes in the fiowgram from 
saturation to divergence with an unstable separatrix in between as a function of coupling 
parameters is a clear sign of TP. It is important to emphasize that the fiowgram method 
does not rely on detecting energy barriers between competing phases. It is simply based on 
phase coexistence typical for first-order transitions. 

We use the following definition of a transition point: for any given set of {g, L) or {V, L) 
we determine the critical value of U from the condition that the ratio of probabilities of 
having zero and non-zero winding numbers is some fixed number C of the order of unity: 

C+ = F(0)/ ^(W-) = C, C^ = F(0)/ = C, (5) 

for the neutral and the charged lines, respectively. We use this particular criterion (among 
many others) because in the thermodynamic limit it gives exact answer for both second- 
and first-order transitions, since the ratio of probabilities in (jSJ is zero in the superfiuid 



phase and infinity in the insulating phase 19|. As mentioned above, at the critical point 



we compute R±. For the scale- invariant continuous transition one expects this quantity to 
approach its universal value as L — ^ oo; for the the I-order transition it has to diverge with 
L since in the superfiuid phase (Wj_) ~ L. The bicritical point is detected by observing 
vanishing charged winding numbers R+ at the neutral line as L is increased, and diverging, 
i?_ oc L, neutral winding numbers on the charged transition line. 

Close to the weakly I-order transition the ratios C± (0) and averages R± exhibit strong 
fiuctuations and are hard to compute precisely. However, fiuctuations in R± are strongly 
correlated with fiuctuations in C±. This property, known as covariance, can be used to 
reduce statistical errors ^Q^. We have employed covariance, by plotting R± vs C±, to reduce 
error bars by nearly an order of magnitude. 

The fiowgram method also allows to "visualize" the renormalization of coupling param- 
eters. If data obtained for {gi^ Li) collapse on top of data for ((72, L2) then one can consider 
this as an evidence that two systems are already in the scaling regime and are governed by 
the same effective coupling at large scales, ges{gi,Li) = 5'ofr(fl'2) -^2)- In practice, we plot 
data for R^ as a function of ln{L/C,{g)) where C,{g) is some interaction-dependent length- 
scale, i.e. we study if system properties at large scales are identical for different coupling 
constants g up to the length-scale renormalization and are part of the same renormalization 
flow. If true, all data points for the DCP model should collapse on a single master curve. 

11 



In the limit of small g, the runaway flow starts as QcsiL) oc ql]^, thus the length scale 
renormalization should follow the C,{g) oc g^^ law. 



Flowgrams for the short-range model 

The most controversial region on the phase diagram for the BH2 model (j21) is located 
in the vicinity of BP, see Fig. |2l where previous simulations failed to resolve the difference 
between BP and upper TP. The standard feature of I-order transitions - a double-peak 
shape of the probability distribution for energy/action ^(5*) - works reliably only if there is 
a signiflcant barrier between two phases. On approach to TP the barrier vanishes and the 
double-peak distribution method fails to work. 

Using flowgrams one can identify the flrst-order transition even if no distinct discontinuity 
in system properties can be clearly seen. The cornerstone of the flrst-order transition is phase 
coexistence. In our case these phases are the insulator characterized by loops of small size 
(in the corresponding channel) and the superfluid in which the number of proliferated loops 
grows with the system size, so that the mean square winding number (W"^) ~ L for large 
enough sizes. The statistical properties of such mixture are drastically different from those of 
the second-order transition point characterized by scaling behavior, where the same number 
of loops of size ~ L must be seen regardless of the (large) system size. The tricritical point 
is then characterized by the unstable separatrix separating the two limiting behaviors. 

In Fig. m we show flows of /?+ vs L for different values of V along the MI-SF boundary. 
The separatrix position ai V = Vltp ~ 0.5 is consistent with lower TP found previously 
lll | . In Fig. the same separatrix is shown to be more pronounced in the neutral channel 
Below we will contrast these flows with the corresponding plots for the DCP action. 

Linear scaling i?+ ~ L is quite obvious close to upper TP as well, where no distinct 
double-peak probability distributions P{S) were seen, see Fig.lHl As the interaction strength 
V crosses the upper tricritical point located at Vjjtp ~ 0.925 ± 0.01 the behavior changes 
from _R+ ~ L to const. 

For the determination of BP, in Fig. [7| we present the flowgram of charged winding 
numbers on the neutral transition line close to the point where the MI-SF transition splits 
into two. For Vltp < < Vbp = 0.90 ± 0.005 we observe diverging behavior typical for 
I-order transitions. For V > Vbp the flow is towards zero values as L increases and is 
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FIG. 4: (Color online). Flowgrams for along the MI-SF line of the short-range model. The 
separatrix is located at lower TP given hy V = Vltp ~ 0.5. 

indicative of the MI-SCF transition where only winding numbers in the neutral, or paired, 
channel proliferate. In contrast, the neutral windings -R_ diverge as ~ L on the charged 
transition line (SCF-SF) because above the BP the SCF superfluid stiffness is finite in the 
SCF phase. This behavior is seen in Fig. |H1 

Fig. El shows neutral winding numbers along the neutral line of MI-SCF transitions 
in the vicinity of BP. The crossover from linear scaling R_ oc L to saturation R_ const 
is very sharp and is associated with passing through BP, where the Il-order MI-SCF line 
terminates at the I-order MI-SF/SCF-SF line. From this plot we estimate Vbp = 0.90±0.005, 
in perfect agreement with the value deduced from flowgrams for the charged line, see Fig. [7| 

Comparing Vbp = 0.90 to Vutp = 0.925 we conclude that upper TP and BP do not 
coincide, and the actual phase diagram (see Fig. ITH)) is topologically identical to the mean 
field result There is, however, a strong (about an order of magnitude) suppression of 
the I-order region above the BP due to strong fluctuation-induced corrections. 
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FIG. 5: (Color online). Flowgrams for R along the MI-SF line of the short-range model. The 
separatrix is located at lower TP given by ^ = Vltp ~ 0.5. 



DECONFINED CRITICALITY 

The field theory of DCP, built on fractionalized spinons (or vortices in the VBS phase), 
hints at the remarkable possibility of a generic continuous SF-VBS transition. Obviously, 
such Il-order transition which occurs between phases characterized by different broken sym- 
metries cannot be derived from the conventional Landau expansion in powers of the cor- 
responding order parameter fields. One cannot even justify the perturbative expansion 
because it should start from the fully disordered quantum groundstate which has no broken 



symmetries (and topological orders). T 
tum groundstate simply does not exist 



lere are strong arguments that disordered quan- 



[2,|l^- Another option is using a phenomenological 
Landau theory which combines superfiuid and solid orders into one multi-component order 
parameter with non-zero modulus as a (non-perturbative) constraint ^ . This approach, 
however, predicts that SF-solid transitions are generically of first-order with the exception 
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FIG. 6: (Color online). Flowgrams for i?+ along the charged SCF-SF line of the short-range model. 
The upper tricritical point is identified from this plot to be at Vutp ~ 0.925. 

of special high-symmetry points. 

At first glance, DCP predictions are in sharp contradiction with the Landau approach. 
However, there is a fundamental difficulty in the DCP theory, namely, the runaway renor- 
malization flow to strong coupling which makes the idea of generic continuous SF-VBS 
transition speculative and without reliable analytical support 7|. In the absence of stable 
perturbative fixed points other possibilities have to be explored as well. One alternative is 
that the DCP action is, in fact, a theory of generic I-order SF-VBS phase transitions 
contrary to the original predictions. Early work on 3D scalar quantum electrodynamics 
did suggest that the renormalization fiow leads to the I-order transition if the number of 
the matter-field components is relatively small (< 100). This prediction, however, is not 
taken seriously because it is known to be incorrect for the one-component system which. 



by the duality mapping 21|, can be shown to have a continuous 3D XY transition. What 
happens in the mult i- component case is a vast open problem in critical phenomena. It 
is worth noting, however, one aspect which is absent in the one-component case and be- 
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FIG. 7: (Color online). Flowgrams for along the neutral MI-SCF line of the short-range model. 
The bicritical point is identified from this plot to be at Vbp = 0.90. 



comes crucial in the case of two and more components — it is the presence of paired phases 
that may intervene between insulator and SF. Pairing effect was shown above to induce 
I-order superfluid-insulator transition for moderate coupling between components. Paring 
correlations in combination with the runaway flow to strong coupling Q would provide a 
natural explanation for the generic weakly I-order SF-VBS transitions in the DCP theory, 
as conjectured in Ref. jl^. 

The phase diagram for the DCP action shown in Fig. |21is not specific to the j— current 
representation. A similar phase diagram is obtained for the XY phase-gauge field formulation 
of the theory Ref. 0, 2^. More refined and extensive simulations have shown that early 



Monte Carlo results for the easy-plane model mentioned in Ref. section IIIA, did not 
go to sufficiently large system sizes to observe deviations from the power-law scaling and 
signatures of the weakly I-order transition Q| . Since for g the length scale for observing 
the I-order transition, if any, is guaranteed to diverge, one has to be extremely careful in 



16 




FIG. 8: (Color online). Flowgrams for i?_ along the charged MI-SF line of the short-range model 
in the SCF phase above the bicritical point (Vbp = 0.90). 

monitoring deviations from the scaling behavior expected for the continuous scale-invariant 
transition. In Refs. P, 0] the second-order SF-VBS transition with the correlation length 
exponent u ^ 0.6 was deduced using straightforward power law fits of the data for system 
sizes L < 24. To make a point of comparison, we studied u for the short-range model (BH2) 
along the SF-MI line for system sizes L < 128 using finite-size scaling of the superfluid 
density derivative, A'^{U = Uc) oc L^/'^~^. We found that (i) the scaling curves are nearly 
perfect straight lines on the log-log plot from the BP all the way to the lower TP, (ii) good 
scaling with v ^ 0.4 was observed up to L = 32 in the middle of the I-order line(!), and 
(iii) in the upper half of the continuous U{1) x f/(l) SF-MI line the slope of the best linear 
fit produced a set of interaction dependent i^iy ) in the range between 0.67 and 0.5. This 
set of data is indicative of the fact that corrections to scaling are extremely important to 
take into account even for continuous transitions. With corrections to scaling included into 
the fits the data become consistent with u ^ 0.67 though the error bars increase by an 
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FIG. 9: (Color online). Flowgrams for along the neutral MI-SCF line of the short-range model. 
This plot is consistent with the existence of BP at Vbp = 0.90 (the same value was deduced from 
flowgrams) and I-order transitions at and slightly below Vbp- 



order of magnitude. We thus conclude that for two-component models v cannot be reliably 
determined on system sizes L < 24. More importantly, one may completely miss the I-order 
transition. 

In what follows we discuss results of worm algorithm Monte Carlo simulations performed 

flQ 

for the DCP action in the self-dual representation |9| (equivalent to the field theory of two 
identical scalar fields with Abelian gauge and U(l)xU(l) symmetry of quartic interactions), 
construct its phase diagram, and prove our earlier assertion that the transition to the SF 
phase is always either I-order or through the SFS state The supersolid state exists 

only for sufficiently large values of coupling to the gauge field g > qbp with gpp ~ 0.6. For 
9 < 9bp we find I-order VBS-SF transition for any finite value of g. 
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FIG. 10: (Color online). Normalized probability density distributions P{S) for the DCP action 
for g = 0.58 and different system sizes; U(L=14)=1.0132, U(L=18)=1.0075, U(L=20)=1.0062, 
U(L=24)=1.0041, U(L=28)=1.00329, U(L=32)=1.00263. 

Bimodal distribution 

Before presenting flowgrams for the DCP action and interpreting them as evidence for 
the I-order deconfinement transition, we would hke to discuss an unambiguous evidence of 
the I-order VBS-SF transition a.t g = 0.58. The significance of this point will become clear 
shortly after we show that this point can be continuously related to smaller values of g by 
the renormalization flow. 

In the vicinity of the VBS-SF critical point deflned according to the criterion (0) we look 
at the probability distribution P{S) and flnetune U to the point where this distribution 
has the most pronounced double-peak structure or flat top. Unfortunately, simulations of 
large L for the long-range action (P) are very time consuming, and U has to be tuned with 
accuracy better then 10~^. Still, we were able to collect reliable statistics for system sizes 
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FIG. 11: Finite-size scaling of the peak positions in Fig. Ilfll 

L < 32, and have observed, see Fig. [TUl the development of the double-peak structure in 
P{S) starting from an anomalously fiat maximum at L = 14. 

We are not aware of any continuous phase transition which features more and more 
pronounced double-peak distribution in energy/action with increasing the system size j^l, 
and consider Fig. to be an unambiguous evidence in favor of the I-order transition. The 
scaling of peak positions with L demonstrating no sign of peaks moving towards each other 
is presented in Fig. ITTl 

Runaway flow and data collapse for small and intermidiate coupling 

In Fig.^Jwe present the DCP action fiowgrams for various coupling constants 0.125 < 
g < 1.2 (the ending point U = 0, (? ~ 1.6436 corresponds to the inverted- XY transition along 
the charged line) and system sizes 8 < L < 24. We observe a remarkable data collapse by 
rescaling system sizes using ^{g) = {g + g"^)^^ for values of (? up to = 0.65 (rescaling is 
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FIG. 12: (Color online). Rescaled flowgrams for i?+ along the charged line. For g < qbp ~ 0.58, 
this line coincides with the VBS-SF deconfinement transition, and, for g > qbp, separates SF and 
SFS phases. System sizes are rescaled hy ^ = {g + g"^)'^ so that the data for g < 0.65 collapse on 
a single master curve. The lines connect data points {g > 0.65) which do not collapse. 



equivalent to the horizontal shift of data points). The flow starts from the universal value 
corresponding to two decoupled U{1) models and, most importantly, reaches the parameter 
range {g = 0.58, L > 14) where we observe the development of the double-peak structure 
in P{S). Several other features of the flow for g < gsp ~ 0.58 are indicative of the I- 
order transition. We observe that Fig. ^1 and Fig. El grow without saturation as 
expected for the mixture of superfiuid and insulating phases. It cannot be reconciled with 
TP below gpp since TP implies an unstable separatrix and precludes data collapse. 

As mentioned in the introduction, data collapse proves that systems with small g and 
linear size L have the same superfiuid response as systems with larger g and smaller L, i.e. 
their long- wave properties are governed by the renormalized coupling gesiOj L). We explicitly 
verify the perturbative renormalization group result ges oc gL of Ref. [l| for small g, and 
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FIG. 13: (Color online) Flowgrams for neutral winding numbers R- along the charged line demon- 
strating growth without saturation at all coupling constants. This behavior is drastically different 
from that observed in Fig. ((S)) for the short-range model, where the separatrix reveals lower TP. 

thus exclude any possibility for having TP at < 0.125. Renormalization procedure allows 
us to extend the flowgram to sufficiently large values of gL > 18 and to relate systems with 
arbitrary small g to the first-order transition because data points for g = 0.58 are collapsing 
on the same continuous curve. [Simulating gL = 18 for g = 0.125 would be equivalent to 
having L = 115 - an impossible task for the long-range model.] Thus, we conclude that the 
I-order transition line goes all the way to the special g = point which is described by the 
f/(l) X f/(l) universality for two independent complex scalar fields 

Flowgrams along the charged line (SFS-SF transition) 

For g > 0.65 the fiow collapse starts to break down. This region features BP and upper 
TP; for g > 0.8 the presence of the paired SFS phase can be seen by the fiowgram method 
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FIG. 14: Flowgrams for charged winding numbers along the charged Une. The same data as 
in Fig. El are now plotted without rescaling to underline linear divergence of R- with system size 
below guTP ~ 0.70. 

even in small samples. The flow collapse on the charged line fails because the order of 
the transition changes at upper TP ( we estimate guTP = 0.70 ± 0.05) which features an 
unstable separatrix. For g > Qutp we are dealing with the continuous SF-SFS transition 
of the inverted- XY type. The diverging flow, then, is replaced with saturation clearly seen 
already for g = 0.75. 

In Fig. El we present data for flowgrams of neutral winding numbers along the same 
(charged) line. This time we do not rescale distances to verify that i?_ obeys the linear 
law oc L characteristic of flnite superfluid response in the SFS and SF phases, or phase 
coexistence at the I-order transition. The divergent flow observed for R_ at all g on the 
charged line by itself is suggesting that BP (where the I-order SF-VBS line splits into SCF- 
SF and SCF-VBS lines) is below guTP ~ 0.70: otherwise, we would observe saturation of 
i?_ between the two milticritical points. Linear scaling of -R+ for g < guTP is presented in 
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FIG. 15: Renormalized flowgrams of R± along the neutral line. For g < gsp ~ 0.58 this line 
coincides with the I-order VBS-SF deconfinement transition and for g > gsp it separates SFS and 
VBS phases. Each symbol corresponds to a value of g shown. The upper set of data points is 
for neutral winding numbers R- (solid symbols), and the lower set is for charged ones i?+ (open 
symbols). The lines connect data points {g > 0.58) which do not collapse on the master curve. 



Flowgrams along the neutral line (VBS-SFS transition) 

Though flowgrams for the charged transition line are sufficient to determine the topology 
of the phase diagram, they are rather insensitive to the location of BP. Instead, one has 
to look at flowgrams for the neutral line to observe how it branches out from the I-order 
SF-VBS line. These are shown in Fig. ^1 A nearly perfect data collapse on divergent flow 
is observed for both types of winding numbers up to (? ~ gsp ~ 0.58 ± 0.02. In Figs. Uni 
and [T7| the same data are plotted without rescaling to underline linear divergence of R± 



Fig.IH 
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FIG. 16: Flowgrams for charged winding numbers along the neutral line. 



with system size below BP. For higher values of g, winding numbers in the charged channel 
start decaying to zero and the flow cannot be collapsed any more. This is an unambiguous 
sign that we have passed BP and the neutral SFS-VBS line has branched out from the I- 
order VBS-SF transition. At the same time, though less dramatically, the neutral winding 
numbers turn to saturation, as expected for the continuous XY transition. 

We summarize our results for the phase diagrams of the DCP and BH2 actions in Fig. 
which is showing fine details invisible in Fig. |2l One may not miss that SF-SFS and SFS- 
VBS lines go extremely close to each other between upper TP and BP. It is only due to 
the flowgram method that we were able to resolve them conclusively. The most important 
result we learn from Fig. is that pairing fluctuations, at the qualitative level correctly 
captured by the mean-field theory, always place the BP below the upper TP regardless of 
the interaction range. 
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FIG. 17: Flowgrams for neutral winding numbers along the neutral line for g > gsp ~ 0.58 
and along the VBS-SF transition for g < gsp -The latter set of the points exhibits linear growth 
due to the I-order transition as in Fig. ()13() . 

CONCLUSION 

To conclude, we have shown that the runaway renormaUzation flow to strong couphng 
in the U(l)xU(l) symmetric two-component scalar quantum electrodynamics in D = 2 + 1 
leads to a I-order deconflnement transition to the SF state. This result, and the structure of 
the phase diagram for relatively strong coupling between the spinous exclude the possibility 
of generic continuous VBS-SF transition in the self-dual deconflned critical action. The 
I-order transition is fluctuation induced — it develops at progressively larger length scales as 
9^0. 

Though in this work we considered a particular DCP model, Eq. (^J, our results have 
a certain degree of universality since the I-order transition was found to be very weak and 
fully developed only in large systems with more than 10*^ sites (all results for small g are 
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FIG. 18: (Color online). Fine structure of phase diagrams for the long-range (left panel), Eq. 
and short-range (right panel), Eq. (j^)), actions in the vicinity of BP. Bold solid lines indicate the 
I-order phase transition. 

universal since they are subject to the renormalization flow). Thus, our study does rule out 
the possibility of having continuous VBS-SF transitions in models with weak-to-intermediate 
long-range coupling. Yet one may argue that our simulations do not rule the possibility of 
continuous VBS-SF transition in some part of the phase diagram of some yet unexplored 
model, say, for large couplings. Of course, they do not. But there is no evidence or even good 
theoretical argument that such a transition line will not be eliminated by the tendency to 
pairing at large g. Such model, if any, would rather be an exception; the generic situation is 
that the combination of the runaway flow and pairing correlations drives a system into either 
the first-order transition for weak and intermediate couplings or to the paired phase for large 
ones. At the phenomenological level, we think that the strongest argument against DCP is 
Landau theory based on the multi-component order parameter with non-zero modulus — so 
far, all numerical data for various microscopic models of the VBS-SF transition are consistent 
with it. 
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